Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology.

BACKGROUND: The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS: A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS: These validated results show significant improvement over current bovine genome annotations.

Generation of a Commercial-Scale Founder Population of Porcine Reproductive and Respiratory Syndrome Virus Resistant Pigs Using CRISPR-Cas.

Disease resistance genes in livestock provide health benefits to animals and opportunities for farmers to meet the growing demand for affordable, high-quality protein. Previously, researchers used gene editing to modify the porcine CD163 gene and demonstrated resistance to a harmful virus that causes porcine reproductive and respiratory syndrome (PRRS). To maximize potential benefits, this disease resistance trait needs to be present in commercially relevant breeding populations for multiplication and distribution of pigs. Toward this goal, a first-of-its-kind, scaled gene editing program was established to introduce a single modified CD163 allele into four genetically diverse, elite porcine lines. This effort produced healthy pigs that resisted PRRS virus infection as determined by macrophage and animal challenges. This founder population will be used for additional disease and trait testing, multiplication, and commercial distribution upon regulatory approval. Applying CRISPR-Cas to eliminate a viral disease represents a major step toward improving animal health.

Ovarian activity and reproductive responses of lactating Angus cows due to a mineral supplementation throughout a timed AI protocol.

Two experiments were conducted to evaluate the effect of intramuscular administration of minerals during a TAI program on the reproductive responses of lactating Angus cows. All cows (n=353) were subjected to a 9-day TAI program based on CIDR insertion plus injections of estradiol, cloprostenol, and eCG, and then TAI 48 h later. In experiment 1, two groups were randomly created, one control with a placebo injection (CON, n=109), and the second received 10 mL of Fosfosan® (MIN, n=172) on day 0 of the synchronization. Conception rate (66.9 vs. 55%) and estrus percentage (55.8 vs. 44%) were higher (P≤0.05) in MIN than in CON cows. Given these results, a second experiment was conducted randomly assigning the cows to two treatments (n=36 each): a single injection of 10 mL of Fosfosan® (MIN-O) on day 0 or two injections of 10 mL of Fosfosan® (MIN-T) on synchronization days 0 and 7. Four cows of each treatment were randomly selected to be scanned with transrectal ultrasound before and during the synchronization protocol to assess ovarian structures and cyclicity, and at day 39 post-TAI for pregnancy diagnosis. Also, blood samples were obtained for the determination of serum minerals and progesterone (P4) concentrations. The number of mineral injections did not affect conception rate (P≥0.1229) conception rate, serum mineral and P4 concentrations, number, and size of emerging follicles, or follicle size according to 1 to 4 classifications. The MIN-T promoted (P<0.05) earlier follicular wave emergence than MIN-O. However, MIN-O cows had a dominant follicle of 15.12 mm, which is more significant (P<0.05) than that in MIN-T cows (13.5 mm). In conclusion, providing a single mineral injection of Fosfosan® at the start of a TAI program is an excellent reproductive strategy in lactating Angus cows to improve the dominant follicle growth, estrus response, and conception rate.

Blood Transcriptomic Analyses Reveal Functional Pathways Associated with Thermotolerance in Pregnant Ewes Exposed to Environmental Heat Stress.

Environmental heat stress triggers a series of compensatory mechanisms in sheep that are dependent on their genetic regulation of thermotolerance. Our objective was to identify genes and regulatory pathways associated with thermotolerance in ewes exposed to heat stress. We performed next-generation RNA sequencing on blood collected from 16 pregnant ewes, which were grouped as tolerant and non-tolerant to heat stress according to a physiological indicator. Additional samples were collected to measure complete blood count. A total of 358 differentially expressed genes were identified after applying selection criteria. Gene expression analysis detected 46 GO terms and 52 KEGG functional pathways. The top-three signaling pathways were p53, RIG-I-like receptor and FoxO, which suggested gene participation in biological processes such as apoptosis, cell signaling and immune response to external stressors. Network analysis revealed ATM, ISG15, IRF7, MDM4, DHX58 and TGFßR1 as over-expressed genes with high regulatory potential. A co-expression network involving the immune-related genes ISG15, IRF7 and DXH58 was detected in lymphocytes and monocytes, which was consistent with hematological findings. In conclusion, transcriptomic analysis revealed a non-viral immune mechanism involving apoptosis, which is induced by external stressors and appears to play an important role in the molecular regulation of heat stress tolerance in ewes.

bESC from cloned embryos do not retain transcriptomic or epigenetic memory from somatic donor cells.

In brief: Epigenetic reprogramming after mammalian somatic cell nuclear transfer is often incomplete, resulting in low efficiency of cloning. However, gene expression and histone modification analysis indicated high similarities in transcriptome and epigenomes of bovine embryonic stem cells from in vitro fertilized and somatic cell nuclear transfer embryos. Abstract: Embryonic stem cells (ESC) indefinitely maintain the pluripotent state of the blastocyst epiblast. Stem cells are invaluable for studying development and lineage commitment, and in livestock, they constitute a useful tool for genomic improvement and in vitro breeding programs. Although these cells have been recently derived from bovine blastocysts, a detailed characterization of their molecular state is lacking. Here, we apply cutting-edge technologies to analyze the transcriptomic and epigenomic landscape of bovine ESC (bESC) obtained from in vitro fertilized (IVF) and somatic cell nuclear transfer (SCNT) embryos. bESC were efficiently derived from SCNT and IVF embryos and expressed pluripotency markers while retaining genome stability. Transcriptome analysis revealed that only 46 genes were differentially expressed between IVF- and SCNT-derived bESC, which did not reflect significant deviation in cellular function. Interrogating histone 3 lysine 4 trimethylation, histone 3 lysine 9 trimethylation, and histone 3 lysine 27 trimethylation with cleavage under targets and tagmentation, we found that the epigenomes of both bESC groups were virtually indistinguishable. Minor epigenetic differences were randomly distributed throughout the genome and were not associated with differentially expressed or developmentally important genes. Finally, the categorization of genomic regions according to their combined histone mark signal demonstrated that all bESC shared the same epigenomic signatures, especially at gene promoters. Overall, we conclude that bESC derived from SCNT and IVF embryos are transcriptomically and epigenetically analogous, allowing for the production of an unlimited source of pluripotent cells from high genetic merit organisms without resorting to transgene-based techniques.

MOSS: multi-omic integration with sparse value decomposition.

SUMMARY: This article presents multi-omic integration with sparse value decomposition (MOSS), a free and open-source R package for integration and feature selection in multiple large omics datasets. This package is computationally efficient and offers biological insight through capabilities, such as cluster analysis and identification of informative omic features. AVAILABILITY AND IMPLEMENTATION: https://CRAN.R-project.org/package=MOSS. SUPPLEMENTARY INFORMATION: Supplementary information can be found at https://github.com/agugonrey/GonzalezReymundez2021.

Host Serum Proteins as Potential Biomarkers of Bovine Tuberculosis Resistance Phenotype.

Eradication of bovine tuberculosis (bTB) continues to be a worldwide challenge. The lack of reliable vaccines dampens the control and eradication programs of Mycobacterium bovis infection and spread. Selection and breeding of cattle resistant to M. bovis infection would greatly enhance the effectiveness of bTB eradication programs. Here, we have evaluated the potential of serum proteins as biomarkers of cattle resistance to bTB in Holstein-Friesian cows, 6-8-year-old, born and raised in similar conditions in herds with bTB prevalence >30%. Serum proteins obtained from uninfected cows (bTB-resistant; R) were compared to those from infected cows (bTB-susceptible; S), defined by a negative or positive bTB diagnosis, respectively. bTB diagnosis included: (i) single intradermal (caudal fold) tuberculin test, (ii) whole blood IFN-gamma test, (iii) gross visible lesions in lymph nodes and lungs by inspection at the abattoir, and (iv) a bacteriological culture for M. bovis. Using 2D-GE and LC-ESI-MS/MS, we found higher expression levels of primary amine oxidase (AO), complement component 5 (C5), and serotransferrin (TF) in R cattle than S cattle. In-house developed and standardized ELISAs for these novel biomarkers showed the best sensitivities of 72, 77, 77%, and specificities of 94, 94, 83%, for AO, C5, and TF, respectively. AUC-ROC (95% CI) values of 0.8935 (0.7906-0.9964), 0.9290 (0.8484-1.010), and 0.8580 (0.7291-0.9869) were obtained at cut-off points of 192.0, 176.5 ng/ml, and 2.1 mg/ml for AO, C5, and TF, respectively. These proteins are involved in inflammatory/immunomodulatory responses to infections and may provide a novel avenue of research to determine the mechanisms of protection against bTB. Overall, our results indicate that these proteins could be novel biomarkers to help identify cattle resistant to bTB, which in turn could be used to strengthen the effectiveness of existing eradication programs against bTB.

Intestinal Transcriptome Analysis Reveals Enrichment of Genes Associated with Immune and Lipid Mechanisms, Favoring Soybean Meal Tolerance in High-Growth Zebrafish (Danio Rerio).

The molecular mechanisms underlying fish tolerance to soybean meal (SBM) remain unclear. Identifying these mechanisms would be beneficial, as this trait favors growth. Two fish replicates from 19 experimental families were fed fishmeal-(100FM) or SBM-based diets supplemented with saponin (50SBM + 2SPN) from juvenile to adult stages. Individuals were selected from families with a genotype-by-environment interaction higher (HG-50SBM + 2SPN, 170 ± 18 mg) or lower (LG-50SBM + 2SPN, 76 ± 10 mg) weight gain on 50SBM + 2SPN for intestinal transcriptomic analysis. A histological evaluation confirmed middle intestinal inflammation in the LG- vs. HG-50SBM + 2SPN group. Enrichment analysis of 665 differentially expressed genes (DEGs) identified pathways associated with immunity and lipid metabolism. Genes linked to intestinal immunity were downregulated in HG fish (mpx, cxcr3.2, cftr, irg1l, itln2, sgk1, nup61l, il22), likely dampening inflammatory responses. Conversely, genes involved in retinol signaling were upregulated (rbp4, stra6, nr2f5), potentially favoring growth by suppressing insulin responses. Genes associated with lipid metabolism were upregulated, including key components of the SREBP (mbtps1, elov5l, elov6l) and cholesterol catabolism (cyp46a1), as well as the downregulation of cyp7a1. These results strongly suggest that transcriptomic changes in lipid metabolism mediate SBM tolerance. Genotypic variations in DEGs may become biomarkers for improving early selection of fish tolerant to SMB or others plant-based diets.

More than one million barriers fragment Europe's rivers.

Rivers support some of Earth's richest biodiversity1 and provide essential ecosystem services to society2, but they are often fragmented by barriers to free flow3. In Europe, attempts to quantify river connectivity have been hampered by the absence of a harmonized barrier database. Here we show that there are at least 1.2 million instream barriers in 36 European countries (with a mean density of 0.74 barriers per kilometre), 68 per cent of which are structures less than two metres in height that are often overlooked. Standardized walkover surveys along 2,715 kilometres of stream length for 147 rivers indicate that existing records underestimate barrier numbers by about 61 per cent. The highest barrier densities occur in the heavily modified rivers of central Europe and the lowest barrier densities occur in the most remote, sparsely populated alpine areas. Across Europe, the main predictors of barrier density are agricultural pressure, density of river-road crossings, extent of surface water and elevation. Relatively unfragmented rivers are still found in the Balkans, the Baltic states and parts of Scandinavia and southern Europe, but these require urgent protection from proposed dam developments. Our findings could inform the implementation of the EU Biodiversity Strategy, which aims to reconnect 25,000 kilometres of Europe's rivers by 2030, but achieving this will require a paradigm shift in river restoration that recognizes the widespread impacts caused by small barriers.

Integrated Transcriptome and Histone Modification Analysis Reveals NDV Infection Under Heat Stress Affects Bursa Development and Proliferation in Susceptible Chicken Line.

Two environmental factors, Newcastle disease and heat stress, are concurrently negatively impacting poultry worldwide and warrant greater attention into developing genetic resistance within chickens. Using two genetically distinct and highly inbred layer lines, Fayoumi and Leghorn, we explored how different genetic backgrounds affect the bursal response to a treatment of simultaneous Newcastle disease virus (NDV) infection at 6 days postinfection (dpi) while under chronic heat stress. The bursa is a primary lymphoid organ within birds and is crucial for the development of B cells. We performed RNA-seq and ChIP-seq targeting histone modifications on bursa tissue. Differential gene expression revealed that Leghorn, compared to Fayoumi, had significant down-regulation in genes involved in cell proliferation, cell cycle, and cell division. Interestingly, we also found greater differences in histone modification levels in response to treatment in Leghorns than Fayoumis, and biological processes enriched in associated target genes of H3K27ac and H3K4me1 were similarly associated with cell cycle and receptor signaling of lymphocytes. Lastly, we found candidate variants between the two genetic lines within exons of differentially expressed genes and regulatory elements with differential histone modification enrichment between the lines, which provides a strong foundation for understanding the effects of genetic variation on NDV resistance under heat stress. This study provides further understanding of the cellular mechanisms affected by NDV infection under heat stress in chicken bursa and identified potential genes and regulatory regions that may be targets for developing genetic resistance within chickens.

T cell epitope content comparison (EpiCC) analysis demonstrates a bivalent PCV2 vaccine has greater T cell epitope overlap with field strains than monovalent PCV2 vaccines.

Porcine circovirus type 2 (PCV2) has one of the highest evolutionary rates among DNA viruses. Traditionally, PCV2 vaccines have been based on the 2a genotype as this was the first genotype discovered. Today, eight genotypes of PCV2 viruses have been identified, and, taken together with the rapid evolutionary rate, propensity to recombine, and high rate of vaccination, further variation in PCV2 is expected. For these reasons, there is a growing genetic gap between available vaccines and field strains. When selecting vaccines, it is important to consider vaccines that contain T cell epitopes that are well-matched to the circulating strains. To quantify the relatedness between PCV2 vaccines and field strains, we predicted and compared their T cell epitope content and calculated Epitope Content Comparison (EpiCC) scores using established in silico tools. T cell epitopes predicted to bind common class I and class II swine leukocyte antigen (SLA) alleles were identified from two major structural proteins, the capsid (encoded by ORF2) and the replicase (encoded by ORF1). The T cell epitope content of three commercial PCV2a-based vaccines (a baculovirus expressed PCV2a ORF2 [VacAlt], a PCV1-PCV2a chimeric virus vaccine [VacA] and a combination cPCV2a-cPCV2b chimeric virus vaccine [VacAB]) and an experimental PCV2b ORF2-based chimeric virus vaccine [VacB] (Table 1), were compared to that of 161 PCV2 field strains (representing genotypes a-f). The T cell epitope content and conservation between vaccine and field strains varied. While all vaccine strains provided broad coverage of the field strains including heterologous genotypes, none of the vaccines covered all the putative T cell epitopes identified in the field strains. PCV2a-based vaccine strains generally scored higher in terms of conserved epitope content against PCV2a field isolates but were not identical. The PCV2b-based vaccine strain had higher scores against PCV2b and PCV2d field strains. The combination PCV2a-PCV2b vaccine (VacAB) had, on average, the highest EpiCC score. PCV2 continues to evolve and EpiCC analysis provides a new tool to assess the possible impact of virus genetic divergence on T cell epitope coverage of vaccine strains. Given that multiple genotypes are currently found and may co-exist on farms, this analysis suggests that a combination of PCV2a and PCV2b vaccine strains may be required to provide optimal coverage of current and future field isolates.

Revisiting the Evolution and Taxonomy of Clostridia, a Phylogenomic Update.

Clostridium is a large genus of obligate anaerobes belonging to the Firmicutes phylum of bacteria, most of which have a Gram-positive cell wall structure. The genus includes significant human and animal pathogens, causative of potentially deadly diseases such as tetanus and botulism. Despite their relevance and many studies suggesting that they are not a monophyletic group, the taxonomy of the group has largely been neglected. Currently, species belonging to the genus are placed in the unnatural order defined as Clostridiales, which includes the class Clostridia. Here, we used genomic data from 779 strains to study the taxonomy and evolution of the group. This analysis allowed us to 1) confirm that the group is composed of more than one genus, 2) detect major differences between pathogens classified as a single species within the group of authentic Clostridium spp. (sensu stricto), 3) identify inconsistencies between taxonomy and toxin evolution that reflect on the pervasive misclassification of strains, and 4) identify differential traits within central metabolism of members of what has been defined earlier and confirmed by us as cluster I. Our analysis shows that the current taxonomic classification of Clostridium species hinders the prediction of functions and traits, suggests a new classification for this fascinating class of bacteria, and highlights the importance of phylogenomics for taxonomic studies.

Identification and quantification of domestic feline cytochrome P450 transcriptome across multiple tissues.

Understanding of cytochrome P450 (CYP) isoform distribution and function in the domestic feline is limited. Only a few studies have defined individual CYP isoforms across metabolically relevant tissues, hampering the ability to predict drug metabolism and potential drug-drug interactions. Using RNA sequencing (RNA-seq), transcriptomes from the 99 Lives Cat Genome Sequencing Initiative databank combined with experimentally acquired whole transcriptome sequencing of healthy, adult male (n = 2) and female (n = 2) domestic felines, expression of 42 CYP isoforms were identified in 20 different tissues. Thirty-seven of these isoforms had not been previously reported in cats. Depending on the tissue, three to twenty-nine CYP isoform transcripts were expressed. The feline genome annotations did not differentiate CYP2E1 and 2E2 genes, demonstrating poor annotation for this gene using the reference genome. As the majority of the sequences are based on automated pipelines, complete cDNA sequences for translation into CYP protein sequences could not be determined. This study is the first to identify and characterize 37 additional CYP isoforms in feline tissues, increasing the number of identified CYP from the previously reported seven isoforms to 42 across 20 tissues.

Polymorphisms within the prolactin and growth hormone/insulin-like growth factor-1 functional pathways associated with fertility traits in Holstein cows raised in a hot-humid climate.

Prolactin (PRL), growth hormone (GH), and insulin-like growth factor-1 (IGF-1) are in hormone-response pathways involved in energy metabolism during thermoregulation processes in cattle. Objective herein was to study the association between single nucleotide polymorphisms (SNP) within genes of the PRL and GH/IGF-1 pathways with fertility traits such as services per conception (SPC) and days open (DO) in Holstein cattle lactating under a hot-humid climate. Ambient temperature and relative humidity were used to calculate the temperature-humidity index (THI) which revealed that the cows were exposed to heat stress conditions from June to November of 2012 in southern Sonora, Mexico. Individual blood samples from all cows were collected, spotted on FTA cards, and used to genotype a 179 tag SNP panel within 44 genes from the PRL and GH/IGF-1 pathways. The associative analyses among SNP genotypes and fertility traits were performed using mixed-effect models. Allele substitution effects were calculated using a regression model that included the genotype term as covariate. Single-SNP association analyses indicated that eight SNP within the genes IGF-1, IGF-1R, IGFBP5, PAPPA1, PMCH, PRLR, SOCS5, and SSTR2 were associated with SPC (P < 0.05), whereas four SNP in the genes GHR, PAPPA2, PRLR, and SOCS4 were associated with DO (P < 0.05). In conclusion, SNP within genes of the PRL and GH/IGF-1 pathways resulted as predictors of reproductive phenotypes in heat-stressed Holstein cows, and these SNP are proposed as candidates for a marker-assisted selection program intended to improve fertility of dairy cattle raised in warm climates.

Neuropeptidome of the Hypothalamus and Pituitary Gland of Indicine × Taurine Heifers: Evidence of Differential Neuropeptide Processing in the Pituitary Gland before and after Puberty.

Puberty in cattle is regulated by an endocrine axis, which includes a complex milieu of neuropeptides in the hypothalamus and pituitary gland. The neuropeptidome of hypothalamic-pituitary gland tissue of pre- (PRE) and postpubertal (POST) Bos indicus-influenced heifers was characterized, followed by quantitative analysis of 51 fertility-related neuropeptides in these tissues. Comparison of peptide abundances with gene expression levels allowed assessment of post-transcriptional peptide processing. On the basis of classical cleavage, 124 mature neuropeptides from 35 precursor proteins were detected in hypothalamus and pituitary gland tissues of three PRE and three POST Brangus heifers. An additional 19 peptides (cerebellins, PEN peptides) previously reported as neuropeptides that did not follow classical cleavage were also identified. In the pre-pubertal hypothalamus, a greater diversity of neuropeptides (25.8%) was identified relative to post-pubertal heifers, while in the pituitary gland, 38.6% more neuropeptides were detected in the post-pubertal heifers. Neuro-tissues of PRE and POST heifers revealed abundance differences ( p < 0.05) in peptides from protein precursors involved in packaging and processing (e.g., the granin family and ProSAAS) or neuron stimulation (PENK, CART, POMC, cerebellins). On their own, the transcriptome data of the precursor genes could not predict the neuropeptide profile in the exact same tissues in several cases. This provides further evidence of the importance of differential processing of the neuropeptide precursors in the pituitary before and after puberty.

Divergent macrophage responses to Mycobacterium bovis among naturally exposed uninfected and infected cattle.

Mycobacterium bovis, the causative agent of bovine tuberculosis (TB), is a successful pathogen that remains an important global threat to livestock. Cattle naturally exposed to M. bovis normally become reactive to the M. bovis-purified protein derivative (tuberculin) skin test; however, some individuals remain negative, suggesting that they may be resistant to infection. To better understand host innate resistance to infection, 26 cattle from herds with a long history of high TB prevalence were included in this study. We investigated the bactericidal activity, the production of reactive oxygen and nitrogen species and the TB-related gene expression profile after in vitro M. bovis challenge of monocyte-derived macrophages from cattle with TB (n=17) and from non-infected, exposed cattle (in-contacts, n=9). The disease status was established based on the tuberculin skin test and blood interferon-gamma test responses, the presence of visible lesions at inspection on abattoirs and the histopathology and culture of M. bovis. Although macrophages from TB-infected cattle enabled M. bovis replication, macrophages from healthy, exposed cattle had twofold lower bacterial loads, overproduced nitric oxide and had lower interleukin (IL)-10 gene expression (P⩽0.05). Higher mRNA expression levels of inducible nitric oxide synthase, C-C motif chemokine ligand 2 and IL-12 were observed in macrophages from all in-contact cattle than in macrophages from their TB-infected counterparts, which expressed more tumour necrosis factor-α; however, the differences were not statistically significant owing to individual variation. These results confirm that macrophage bactericidal responses have a crucial role in innate resistance to M. bovis infection in cattle.

Identification of Gene Networks for Residual Feed Intake in Angus Cattle Using Genomic Prediction and RNA-seq.

Improvement in feed conversion efficiency can improve the sustainability of beef cattle production, but genomic selection for feed efficiency affects many underlying molecular networks and physiological traits. This study describes the differences between steer progeny of two influential Angus bulls with divergent genomic predictions for residual feed intake (RFI). Eight steer progeny of each sire were phenotyped for growth and feed intake from 8 mo. of age (average BW 254 kg, with a mean difference between sire groups of 4.8 kg) until slaughter at 14-16 mo. of age (average BW 534 kg, sire group difference of 28.8 kg). Terminal samples from pituitary gland, skeletal muscle, liver, adipose, and duodenum were collected from each steer for transcriptome sequencing. Gene expression networks were derived using partial correlation and information theory (PCIT), including differentially expressed (DE) genes, tissue specific (TS) genes, transcription factors (TF), and genes associated with RFI from a genome-wide association study (GWAS). Relative to progeny of the high RFI sire, progeny of the low RFI sire had -0.56 kg/d finishing period RFI (P = 0.05), -1.08 finishing period feed conversion ratio (P = 0.01), +3.3 kg^0.75 finishing period metabolic mid-weight (MMW; P = 0.04), +28.8 kg final body weight (P = 0.01), -12.9 feed bunk visits per day (P = 0.02) with +0.60 min/visit duration (P = 0.01), and +0.0045 carcass specific gravity (weight in air/weight in air-weight in water, a predictor of carcass fat content; P = 0.03). RNA-seq identified 633 DE genes between sire groups among 17,016 expressed genes. PCIT analysis identified >115,000 significant co-expression correlations between genes and 25 TF hubs, i.e. controllers of clusters of DE, TS, and GWAS SNP genes. Pathway analysis suggests low RFI bull progeny possess heightened gut inflammation and reduced fat deposition. This multi-omics analysis shows how differences in RFI genomic breeding values can impact other traits and gene co-expression networks.

A clone-free, single molecule map of the domestic cow (Bos taurus) genome.

BACKGROUND: The cattle (Bos taurus) genome was originally selected for sequencing due to its economic importance and unique biology as a model organism for understanding other ruminants, or mammals. Currently, there are two cattle genome sequence assemblies (UMD3.1 and Btau4.6) from groups using dissimilar assembly algorithms, which were complemented by genetic and physical map resources. However, past comparisons between these assemblies revealed substantial differences. Consequently, such discordances have engendered ambiguities when using reference sequence data, impacting genomic studies in cattle and motivating construction of a new optical map resource--BtOM1.0--to guide comparisons and improvements to the current sequence builds. Accordingly, our comprehensive comparisons of BtOM1.0 against the UMD3.1 and Btau4.6 sequence builds tabulate large-to-immediate scale discordances requiring mediation. RESULTS: The optical map, BtOM1.0, spanning the B. taurus genome (Hereford breed, L1 Dominette 01449) was assembled from an optical map dataset consisting of 2,973,315 (439 X; raw dataset size before assembly) single molecule optical maps (Rmaps; 1 Rmap = 1 restriction mapped DNA molecule) generated by the Optical Mapping System. The BamHI map spans 2,575.30 Mb and comprises 78 optical contigs assembled by a combination of iterative (using the reference sequence: UMD3.1) and de novo assembly techniques. BtOM1.0 is a high-resolution physical map featuring an average restriction fragment size of 8.91 Kb. Comparisons of BtOM1.0 vs. UMD3.1, or Btau4.6, revealed that Btau4.6 presented far more discordances (7,463) vs. UMD3.1 (4,754). Overall, we found that Btau4.6 presented almost double the number of discordances than UMD3.1 across most of the 6 categories of sequence vs. map discrepancies, which are: COMPLEX (misassembly), DELs (extraneous sequences), INSs (missing sequences), ITs (Inverted/Translocated sequences), ECs (extra restriction cuts) and MCs (missing restriction cuts). CONCLUSION: Alignments of UMD3.1 and Btau4.6 to BtOM1.0 reveal discordances commensurate with previous reports, and affirm the NCBI's current designation of UMD3.1 sequence assembly as the "reference assembly" and the Btau4.6 as the "alternate assembly." The cattle genome optical map, BtOM1.0, when used as a comprehensive and largely independent guide, will greatly assist improvements to existing sequence builds, and later serve as an accurate physical scaffold for studies concerning the comparative genomics of cattle breeds.

Transcriptional Response to Acute Thermal Exposure in Juvenile Chinook Salmon Determined by RNAseq.

Thermal exposure is a serious and growing challenge facing fish species worldwide. Chinook salmon (Oncorhynchus tshawytscha) living in the southern portion of their native range are particularly likely to encounter warmer water due to a confluence of factors. River alterations have increased the likelihood that juveniles will be exposed to warm water temperatures during their freshwater life stage, which can negatively impact survival, growth, and development and pose a threat to dwindling salmon populations. To better understand how acute thermal exposure affects the biology of salmon, we performed a transcriptional analysis of gill tissue from Chinook salmon juveniles reared at 12° and exposed acutely to water temperatures ranging from ideal to potentially lethal (12° to 25°). Reverse-transcribed RNA libraries were sequenced on the Illumina HiSeq2000 platform and a de novo reference transcriptome was created. Differentially expressed transcripts were annotated using Blast2GO and relevant gene clusters were identified. In addition to a high degree of downregulation of a wide range of genes, we found upregulation of genes involved in protein folding/rescue, protein degradation, cell death, oxidative stress, metabolism, inflammation/immunity, transcription/translation, ion transport, cell cycle/growth, cell signaling, cellular trafficking, and structure/cytoskeleton. These results demonstrate the complex multi-modal cellular response to thermal stress in juvenile salmon.

Expression of autocrine prolactin and the short isoform of prolactin receptor are associated with inflammatory response and apoptosis in monocytes stimulated with Mycobacterium bovis proteins.

Increased levels of prolactin (PRL) have recently been associated with carcinogenesis and the exacerbation of autoimmune diseases, and might be involved in the progression of tuberculosis (TB). To investigate the relationship between PRL and prolactin receptor (PRLr) expression with inflammatory response and apoptosis in monocytes, we used THP-1 cells stimulated with antigens of the Mycobacterium bovis AN5 strain culture filtrate protein (CFP-M. bovis). Western blot (WB), real-time Polymerase chain reaction (PCR), and immunocytochemistry were performed to identify both PRL and PRLr molecules. PRL bioactivity and proinflammatory cytokine detection were assessed. The results showed that PRL and PRLr messenger RNA (mRNA) were synthesized in THP-1 monocytes induced with CFP-M. bovis at peaks of 176- and 404-fold, respectively. PRL forms of 60 and 80kDa and PRLr isoforms of 40, 50, and 65kDa were also identified as time-dependent, while 60-kDa PRL, as well as 40-, and 50-kDa PRLr, were found as soluble forms in culture media and later in the nucleus of THP-1 monocytes. PRL of 60kDa released by monocytes exhibited bioactivity in Nb2 cells, and both synthesized PRL and synthesized PRLr were related with nitrite and proinflammatory cytokine levels proapoptotic activity in CFP-M. bovis-induced monocytes. Our results suggest the overexpression of a full-autocrine loop of PRL and PRLr in monocytes that enhances the inflammatory response and apoptosis after priming with M. bovis antigens.